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Abstract. The space L of oriented lines, or rays, in is a 4-dimensional 
space with an abundance of natural geometric structure. In particular, it 
boasts a neutral Kahler metric which is closely related to the Euclidean metric 
on . In this paper we explore the relationship between the focal set of a line 
congruence (or 2-parameter family of oriented lines in M.^) and the geometry 
induced on the associated surface in L. The physical context of such sets is 
geometric optics in a homogeneous isotropic medium, and so, to illustrate the 
method, we compute the focal set of the fc*"^ reflection of a point source off the 
inside of a cylinder. The focal sets, which we explicitly parameterize, exhibit 
unexpected symmetries, and are found to fit well with observable phenomena. 



The space of oriented afSne lines in , which we denote L, has an abundance of 
natural geometric structure. The purpose of this paper is to continue recent work 
[7] relating this structure to geometric optics in a homogeneous isotropic medium: 
the theory of light propagation under the assumption that the light travels along 
straight lines in R^. The fundamental objects of study are 2-parameter families of 
oriented lines, or line congruences^ which we view as surfaces in L. Thus we are 
lead to consider the geometry of immersed surfaces S C L. 

In the first instance, since L can be identified with the tangent space to the 
2-sphere, there is the natural bundle map tt : L — > 5^. If ttIs : S ^ is not an 
immersion, we say that S is flat. Otherwise, E can be described, at least locally, 
by sections of the canonical bundle. 

On the other hand, there is a natural symplectic structure on L, and S C 
L is lagrangian with respect to this symplectic structure iff the line congruence 
admits a family of orthogonal surfaces in M.^. In geometric optics such surfaces 
are the wavefronts of the propagating light. As a wavefront evolves along the line 
congruence, if there is any focusing, the surface becomes singular. The points at 
which this occurs are referred to as focal points. 

In addition, L admits a natural complex structure JJ, which, together with the 
symplectic structure, forms a natural Kahler structure [8] . The metric G is of signa- 
ture (+H ) and therefore the metric induced on a surface E may be riemannian, 

lorentz or degenerate. 

The aim of this paper is two-fold: to relate the geometry induced on S C L by 
Gr with the set of focal points of the line congruence in R'^, and to demonstrate the 
effectiveness of the geometric formalism in computing focal sets explicitly. More- 
over, we hope to demonstrate that the results of such computations can have visible 
physical significance. 

For the first aim we prove: 
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Main Theorem 1. 

Let S be an immersed surface in L. // S is flat, there is exactly one focal point 
on each line of the congruence. If S is not flat, then there is none, one or two focal 
points on each line iff the metric induced on 12 by G is riemannian, degenerate or 
lorentz (respectively). 

In addition, we relate the distance between pairs of focal surfaces and the angle 

between their normals to geometric quantities on E C L. 

Mathematically, the focal set of a generic line congruence is well understood [1] 
[3] [11]. Special examples of focal sets, also referred to as caustics by some authors, 
have been studied for many decades [2] [4] [5] [10]. For the second aim, we compute 
explicitly the focal set of the line congruence formed by the multiple reflection of a 
point source off the inside of a cylinder. The first such reflected focal set is often 
referred to as the coffeecup caustic, since its cross-section is commonly observed on 
the top of a cup of coffee in the presence of a bright light. We prove: 

Main Theorem 2. 

Consider the reflection of a point source at (—1, 0,0) off the inside of a cylin- 
der lying along the a? -axis with radius a. The focal set of the reflected line congru- 
ence consists of a surface: 

— ( 1)^^+^/ [/sini> i + (a^ — P sin^ v)^'^''[2klcosvsmv e~™ + (a^ — P sin^ v)^\ 

a?'^ [2fc/ cos u + (a^ — P sir? u) 2 ] 

3 _ k{l-u'^)[a'^ -t^ -2l'^s\x?v + 2klcosv{a'^ -l'^s\r?v)^ 
u[2A;Zcost; + (a^ — P siv? v)^] 
and a curve in the x^x^ -plane: 

z = {-lf+^ka-^^[lsmv i + {a^- l"^ six? vY^f^[l + 2ke'^a^ - l^sii?v)^ 
where z = x^ -\- io? , u e K and v is in the domain 

< V < IT for I < a 

and 

— sm~^{a/l) < V < sm~^{a/l) for I > a. 



Thus the focal set consists of two sets: a translation invariant surface and a curve 
lying in the plane perpendicular to the symmetry axis and containing the source 
point. The commonly observed coffeecup caustic is the level sets of the former with 
A: = 1, while the latter lies outside the cup and so is not observed. With the aid 
of a polished brass cylinder, the first author, in collaboration with Grace Weir, has 
photographed multiple reflections of up to 4*^ order. The results agree closely with 
the result of Main Theorem 2. 

This paper is organised as follows. The next section contains a review of the 
complex geometric structure on the space of oriented lines in R^, as developed 
in [8]. The second section describes the first order invariants of an arbitrary line 
congruence, while the next section contains an exposition on the construction of 
surfaces using line congruences. This is required in order to have a description 
of the reflective surfaces. To illustrate the method we compute the normal line 
congruence to the parabaloids in M^. 
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Section 4 turns to reflection in a surface, as expounded more fully in [7]. By 
way of example, we compute the reflection of a plane wave off the parabaloids 
computed in the previous section. In section 5 we describe the general procedure 
for computing the focal set of an arbitrary line congruence and prove Main Theorem 
1. 

Finally we look in detail at the coffeecup caustics: the focal sets formed by 
reflection of a point source off the inside of a cylinder. We do this for multiple 
reflections off the cylinder and find that the focal sets consists of two distinct sets: 
one a translation invariant surface and the other a curve outside of the cylinder (and 
hence not physically seen). In the limit as the source tends to infinity, the incoming 
wave is flat and the resulting focal set is independent of the angle of incidence of 
the wave. Section 7 contains a discussion of these focal sets and their properties. 

1. The Space of Oriented lines 

We start with 3-dimensional Euclidean space M'' and flx standard coordinates 
{x^,x'^,x^). In what follows we combine the first two coordinates to form a single 
complex coordinate z = + ix'^, set t = x^ and refer to coordinates (z, t) on R'^. 

Let L be the set of oriented lines, or rays, in Euclidean space R'^. Such a line 7 
is uniquely determined by its unit direction vector U and the vector V joining the 
origin to the point on the line that lies closest to the origin. That is, 

7 = {V + rlJeR^|reR} 

where r is an affine parameter along the line. 

By parallel translation, we move U to the origin and V to the head of U. Thus, we 
obtain a vector that is tangent to the unit 2-dimensional sphere in R^ . The mapping 
is one-to-one and so it identifies the space of oriented lines with the tangent bundle 
of the 2-sphere T S^ (see Figure 1). 

L = { (U,V) e R^ X R=^ I |IJ| = 1 IJ-V = 0} 
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FIGURE 1 
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L is a 4-dimensional manifold and the above identification gives a natural set 
of local complex coordinates. Let ^ be the local complex coordinate on the unit 
2-sphere in obtained by stereographic projection from the south pole. 

In terms of the standard spherical polar angles {0,(j)), we have ^ = tan(|)e"'^. 
We convert from coordinates ^) back to {6, (j)) using 



cos6i=i^ sin6'=^^ cos(?!)=^±i sin<?!)=^^ 



(1.1) 



This can be extended to complex coordinates r]) on L minus the tangent space 

over the south pole, as follows. First note that a tangent vector X to the 2-sphcrc 
can always be expressed as a linear combination of the tangent vectors generated 
by 9 and (f): 

In our complex formalism, we have the natural complex tangent vector 
^ cos^(f)^^ ' 



5C \d9 2cos(f)sin(f)a<^y 

and any real tangent vector can be written as 

^ d d 

for a complex number 77. We identify the real tangent vector X on the 2-sphere 
(and hence the ray in M?) with the two complex numbers (^, rj). Loosely speaking, 
^ determines the direction of the ray. and 77 determines its perpendicular distance 
vector to the origin - complex representations of U and V. 

The coordinates (^, rj) do not cover all of L - they omit all of the lines pointing 
directly downwards. However, the construction can also be carried out using stereo- 
graphic projection from the north pole, yielding a coordinate system that covers all 
of L except for the lines pointing directly upwards. Between these two coordinate 
patches the whole of the space of oriented lines is covered. In what follows we work 
in the patch that omits the south direction. 

While the coordinates depend on a choice of origin, a translation of the origin 
simply alters the coordinates by : 

^^C' = C ri^r,' = ri + a + h^-ae, (1-2) 
while a rotation about the origin is given by 

for a, /3 e C and a, 6 e IR with a? +(3(3 = 1. 

Definition 1. The map $ : L x M ^ is defined to take ((^,?7),r) G L x R to 
the point in on the oriented line (^, rj) that lies a distance r from the point on 
the line closest to the origin (see the Figure 2). 
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This map, which is of crucial importance when describing surfaces in R^, has the 
following cooordinate expression: 

Proposition 1. [6] //$(^,?7,r) = (z(^, 77, r), t(^, r;, r)), then: 



2(r?-77g2)+2j(l + eg)r 



t = 



(1.4) 



where z ^ + ix^ , t = and (x^ , x^ , x^ ) are Euclidean coordinates in MP. 

The existence of complex coordinates implies that there is a complex structure 
J on L. In fact, this map J : T^L —f T^L can be defined as follows: given a line 7 
in R'^, the tangent space T^L can be identified with the Jacobi fields along the line 
that are orthogonal to the direction of 7. Then J is given by rotation of the Jacobi 
field through 90*^ about the line 7 [9] . 

The complex structure can be supplemented with a natural symplectic structure 
ri. This is obtained by pulling back the canonical symplectic structure on T*S^ with 
the round metric, considered as a mapping from T to T*S^. This symplectic 
structure is compatible with the complex structure and so together they define a 
metric by G(-,-) = n{I-,-). 

The metric G has many interesting properties: it is of signature (+4 ), con- 
formally flat and scalar flat (but not Einstein). In addition, the identity component 
of the isometry group of G is isomorphic to the identity component of the Euclidean 
isometry group [8]. A Kahler structure with this property on the space of oriented 
lines in R" exists only when n = 3 or n = 7 [13]. The metric has the following inter- 
pretation: the norm of a vector in T^L with repsect to G is the angular momentum 
of the associated Jacobi field along 7. 

In the next section we investigate the geometric structures induced by (J,f7,G) 
on a surface in L. 
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2. Line Congruences 

In its simplest form, geometric optics models the propagation of light through a 
homogeneous isotropic medium by a 2-parameter family of rays. 

Definition 2. A line congruence is a 2-parameter family of oriented lines in M^. 

Prom our perspective a line congruence is a surface S in L. For example, a 
point source corresponds to the 2-paramcter family of oriented lines that contain 
the source point, which thus defines a 2-sphcrc in L. 

For computational purposes, we must give explicit local parameterizations of 
the line congruence. In practice, it will be given locally by a map C ^ L : i— > 
{^{n,p,),r]{n,p,)). A convenient choice of parameterization will often depend upon 
the specifics of the situation, but our formalism holds for arbitrary parameteriza- 
tions. 

The dual picture of light propagation is to consider the wavefronts, or surfaces 
that are orthogonal to a given set of rays. However, not every line congruence has 
such orthogonal surfaces - indeed, most don't. To explain this wo consider the first 
order properties of S, which can described by two complex functions, the optical 
scalars: p, cr : S x M — > C The real part 6 and the imaginary part A of p are the 
divergence and twist of the congruence, while a is the shear. 

Definition 3. A null frame in is a trio {eo, e+, e_} of complex vector fields in 
C (?) TK"^, where eo is real, e+ is the complex conjugate of e_ and they satisiy the 
following orthogonality properties: 

<eo,eo>=l <eo,e+>=0 <e+,e+>=0 <e+,e_>=l, 

where we have extended the Euclidean inner product < •, • > of bilinearly over 
C. Orthonormal frames {cq, ei, 62} on TM'^ and null frames are related by 

e+ = -^(ei -iea) e_ = -^(ei ies). (2.1) 

Definition 4. A congruence null frame for S c L is a null frame {eo,e+,e_} if, 
for each 7 € S, we have eo tangent to 7 in M^, and the orientation of {eo, ei, 62} is 
the standard orientation on M^. 

Definition 5. Given a line congruence and null frame, the optical scalars are 
defined by: 

P =< Veo 6+, e_ > (J =< Veo e+, e+ > 

where V is the Euclidean connection on ffi.^. 

These have the following geometric interpretation. Consider a specific ray 7 in 
the line congruence and a point p along this ray. Now consider the unit circle in the 
plane orthogonal to the ray at p. As we flow this circle along the line congruoncc 
this circle will become distorted. To first order in the affine parameter along the ray 
the real part of p measures the divergence or contraction of the circle, the imaginary 
part determines the rotation of the circle, while a measures the shearing [12] (see 
Figure 3). 
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FIGURE 3 



For computational purposes the optical scalars are given by: 

Proposition 2. [6] For a parameterized line congruence the optical scalars have 
the following expressions in terms of first derivatives of the parameterization: 



d^rid i — d rid^ — d~r] d £ 

p = + \i = 



(2.2) 



d rjd T] — d^rjd^rj d rjd rj — d^rjd^rj 

where 

o ri = an + rat, = a n = ori + rot =, 

and d and d are differentiation with respect to fi and p,, respectively. 

The twist has the following important interpretation: 

Proposition 3. There exists surfaces orthogonal to the rays of a line congruence 
if and only if the twist of the line congruence is zero. 

Moreover, in terms of the Kahler structure, the optical scalars have the following 
significance: 

Theorem 1. [8] 

A line congruence E C L is lagrangian (i.e. fll^: — 0) iff the twist of Yi is zero. 

A line congruence S C L is holomorphic (i.e. J preserves the tangent space TYi) 
iff the shear of E is zero. 

The metric induced on S by G is riemannian, degenerate or lorentz iff \o^\'^ < A^, 
I dp = or |(Tp > \^ , respectively, where A is the imaginary part of p. 

A further geometric quantity is the curvature of the line congruence, which is 
defined to he k = pp — aa. A line congruence will be said to be flat if k = 0. 
If a line congruence is non-flat, then the direction of the congruence can be used 
as a parameterization [6]. In other words the line congruence is locally given by 
C The point source line congruence is non-flat, while the set of rays 

orthogonal to a given line forms a flat line congruence. 
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3. Constructing Surfaces Using Line Congruences 

We now describe how to construct surfaces in R"^ using line congruences. Given 
a line congruence E C L, a map r : E — > R determines a map E — > R"^ by 77) i— )■ 
<i>((^, ry), r(^, 77)) for (^,77) G E. In other words, we pick out one point on each line 
in the congruence (see Figure 4). 

With a local parameterization /i of E, composition with the above map yields a 
map C — + R"^ which comes from substituting r = r{fi, p.) in equations (1.4). 



J 






1 \\r=r(|x,|x) 


FIGURE 4 



Of particular interest are the surfaces in M'^ orthogonal to the line congruence - 
when the line congruence is normal. As mentioned earlier, these exist iff the twist 
of the congruence vanishes. By the first of equation (2.2), this is an integrability 
condition for a real function: 

Theorem 2. [6] 

A line congruence {^{^, fl),ri{^, p.)) is orthogonal to a surface in R"^ iff there 
exists a real function r{fi, p,) satisfying: 

If there exists one solution, there exists a 1-parameter family generated by a real 
constant of integration. An explicit parameterization of these surfaces m R^ is given 
by inserting r — r(/i,/i) in (1.4)- 

We now show how to construct the normal line congruence of a non-flat oriented 
surface S" C R'^. As it is non-flat, we can parameterize S by the direction ^ of 
the unit normal. Equivalently, we use the real stereographic projection coordinates 
(0,0). 

Suppose S is defined by the pre-image of zero of a function G : R^ ^ R. Then 
the unit normal to S is, using standard Euclidean coordinates {x^ , x'^ , x'^) , 
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where 

(\ -1/2 
^/5G\M dG 

The Unk with the coordinates ^ = tan(|)e'"''' is 

Ai(a;^, a;^, a;^) = cos<?f)sin^ ^2(3^^, a;^, a;^) = sin^sin^ (3.2) 

Az{x^,x'^,x^) ^ cose (3.3) 

The tangent plane through a point ((?, (j)) on S is given by 

Aix^ + A^x^ + Agx^ = B (3.4) 

where Ai are given by (3.2) to (3.3) and B is a function of (0, 0) determined by the 
surface. 

The function r(^,^) is the distance of the point on the surface to the point on 
the normal line which lies closest to the origin. This is given by 



^Al + Al + Ai 

Finally, the exact functional relationship between ^ and r] is given by 

— 1 — (rV 

The task then, reduces to finding r as a function of 9 and (f), or equivalently, 

^ and ^. In many simple cases it is possible to solve (3.2) to (3.3) and express 
{x^ ,x^, x^) as functions of 9 and (j). These can be directly inserted into (3.4) to find 
r{9A)- 

Once we have r and 77 as functions of ^, equations (1.4) give the explicit param- 
eterization of S in terms of ^. 

We now work through an example in detail. 

Example: Elliptic and Hyperbolic Paraboloids 
Suppose the surface S is determined by 

a 

for some constants a and h. The elliptic paraboloid with a = 1 6 = 1 and the 
hyperbolic paraboloid with a = —b = 1 are graphed below. 
The unit normal is 

\ b'^ J \ a dx^ b dx^ dx^ j 

We can invert the the relations (3.2) to (3.3) to 

■I a , . lb , _ o f 0- 1 I b , o \ 2^ 
— cos (b tan 9 x = - cos (b tan 9 x = — \ - cos + - sm 6 \ tan 9 



Note that the coordinate domain O<0<7r/2,O<^<27r cover all of the 
paraboloid. 
These give 



B = — {a cos^ (j) + b siv? 



sin^ I 



4cos^ 
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Finally converting to holomorphic coordinates to get 

' 4(1-^0(1+^0 4(1-^0' 

The equations (1.4) now give the explicit parameterization of the paraboloid. 
Figure 5 shows the resulting parameterization of the parabaloids with a = —b = 1 
and a ^ b = I. 




The lines on these surface are the images of the lines of longitude and latitude 
around the north pole under the coordinates. 

4. Reflection 

We turn now to the reflection of an oriented line in a surface in MP. This 
is equivalent to the action of a certain group on the space of oriented lines, as 
described by: 

Theorem 3. [7] 

Consider a parametric line congruence ^ — ^i(/ii,/ii), rj = ?7i(/ii,/2i) reflected 
off an oriented surface with parameterized normal line congruence ^ — £,o{^J.Q, p-o), 
V = %(Mo,A''o) ond r = 'ro(/ioj/^o) satisfying (3.1) with ^ = c*"''^ V = Vo- Then 
the reflected line congruence (^2,V2) is 

(i-eoeoKi-2eo' 

„o = (^0-6)" _ (l+go6)' ^ , (|"o-6)(l+Co6)(l+$ol"o) „ 

'2 ((l-?oCo)?i-2?o)^ '0 ((l-?o«o)«i-2«o)^ ((l-?o?o)?i-2?o)2 '0' ^ ^ 

(4.2) 

where the incoming rays are only reflected if they satisfy the intersection equation 
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By virtue of the intersection equation, an alternative way of writing (4.2) is 

„ _ -{l+ioiof , 2(4~o-6)(l+Co6)(l+CoCo) 

The geometric content of this is: reflection of an oriented line can be decomposed 
into a sum of a translation and a rotation about the origin. 

Example: Plane wave reflected off a paraboloid 

Consider the paraboloid, as given earlier in section 3. Thus, it is parameterized 
by the direction of the normal, with 

4(1-^0^0)2 

and 

a(5o+&)'-&(eo-eo)' 



To = 



4(l-?o^o)(l+&eo) 



Assume that the incoming plane wave has normal direction along the positive 
a;^-axis, that is, = 1. By the reflection equation (4.1), the resulting direction is 

2^0 + 1 - eoCo 
l-^o^o-2fo' 

and substituting the equation of the paraboloid into (4.2) yields 



(l-Co6)2(l-eo6- 2^0)2 
K^o - ^o)(2 + 36 - 3^0 + 2eg + 2gg - 2goCo + 2Co^g - 2Cggo) 
(1 - 6^0)2(1 -eo^o- 2^0)2 
^'(^0 - ?"o)(2Co'C"o - 2^0^! - 2Co'eo - S^o'el + ^^Mo " ^^Mo) 



(4.5) 



+ 



(1 - 6^0)2(1 - Cofo - 2Co)2 



Since the incoming ray direction is fixed, we are parameterizing the reflected 

line congruence by the direction ^ of the normal to the surface at the point of 
reflection. A direct integration of equation (3.1) gives the function r as 

2(6 + Co) [a(eo + Co)^ - bj^o - 6)2] ^ 
(1 + 66)2(1-66) 

Finally, the wavefronts form a one-parameter family of parameterized surfaces, 
which can be obtained by substituting for 6 V and r in equation (1.4). The result 
in spherical polar coordinates, after some simplifications, boils down to: 

= 8a(l - cos e) tan 6 cos^ + C(l - 2 sin^ 6 cos^ (p) 
= 2(2a sin^ 9 cos^ + 6) tan 9 sin (f> — 2C sin^ sin (j) cos (f> 

= a(4 cos^ 9-1) tan^ 9 cos^ </> - 6 tan^ 6* sin^ - 2C sin 6* cos 9 cos ^ 
where < 6* < 7r/2 and < </)< 27r. 
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Figure 6 shows examples of the reflected wavefronts for an elliptic and a hyper- 
bolic paraboloid. 

5. The Focal Set of a Line Congruence 
Let a and p be the optical scalars of a line congruence S as described in section 

2. 

Definition 6. A point p on a line 7 in a line congruence is a jocal 'point if p and 
(J blow-up at p. The set of focal points of a line congruence S generically form 
surfaces in R"^, which will be referred to as the jocal surfaces of S. 

Theorem 4. The focal set of a parametric line congruence S is 

{ $(7, r) I 7 e S and 1 - (po + Po)r + (poPo - o-octo)?"^ = 0}, 

where the coefficients of the quadratic equation are given locally by (2.2) evaluated 
atr = 0. 

Proof. In terms of the afhne parameter r along a given line, the Sachs equations, 
which a and p must satisfy, are [12]: 

dp 2 — da 

These are equivalent to the vanishing of certain components of the Ricci tensor of 
the Euclidean metric. They have solution: 

_ Po - jpoPo - cro'ao)r ^ _ oq 

1 - (po + po)r + (poPo - croCTo)r2 1 - (po + po)r + (poPo - (Joao)r^ ' 

where cfq and po are the values of the optical scalars at r = 0. The theorem 
follows. □ 

This has the following corollary: 

Corollary 1. Let Yi he a line congruence, p — Q -\- Xi, a the associated optical 
scalars and po, Oo, Ao, co their values at r = 0. 
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// S is flat with non-zero divergence, then there exists a unique foeal surface S 
given by r — {2Qq)^^ . If it is fl,at with zero divergence, then the focal set is em,pty. 

IfTi is non-flat, then there exists a unique focal point on each line iff |croP = Aq, 
there exist two focal points on each line iff\cTo\^ < Aq and there are no focal points 
on each line iff\cFo\'^ > Aq. The focal set is given by 

^_ eo±(Kp-Ag)i _ 

PoPo — crocro 

Proof. The focal set of a parameterized line congruence are given by r = r{iJ,,p,) 
satisfying the quadratic equation in Theorem 4. If k = 0, then there is none or one 
solution depending on whether Oq = or not. 

If K ^ then there are two, one or no solutions iff |crop — Aq is greater than, 
equal to or less than zero (respectively). 

The solution of the quadratic equation in each case is as stated. □ 

There is also the equivalent definition for focal surfaces: 

Proposition 4. A continuously differentiable surface S in is a focal surface of 
a line congruence E iff every line in E is tangent to S at some point. 

Proof. Let S be locally parameterized by /x (C(a*; m)) ^(a*) m))) where t]) are the 
canonical coordinates above. Then a surface in given by r : S — > M is tangent 
to the line congruence iff 



Det 



2C 2C i-;C 
i+«f i+(J i+« 
dz dz dt 



0, 



dz dz dt 

where the partial derivatives are in fj, and Jl. This determinant equation, which is 
2^{dzdt - dzdt) + 2^{dzdt - dzdt) + (1 - £,C){dzdz - dzdz) = 0, 

is a quadratic equation for r = r(/i, jl) with coefficients given by the first derivatives 

of z{^,fi), z{^,fi) and t(p,fi). 

Carrying out the differentiation we find, for example, that 

a_ _ / 2r mv+in) \ai: , / 2(l+2g|)CV _ Mv-efi) \aE 



' (i+«)' ' ' i+«r 

Similar computations finally yield the quadratic that appears in Theorem 4. □ 

We now explore the geometric properties of the focal set. First we prove: 
Main Theorem 1. 

Let S be an immersed surface in L. // S is flat, there is exactly one focal point 
on each line of the congruence. IfHis not flat then there is none, one or two focal 
points on each line iff the metric induced on S by G is riemannian, degenerate or 

lorentz (respectively). 

Proof. The number of focal points on a given line is determined by the sign of the 
discriminant of the quadratic equation in Theorem 4: |fTo|^ — Ag. By Theorem 1 (c/. 
Theorem 2 of [8]), this is precisely what determines the sign of the metric induced 
on S by G: the metric is riemannian, degenerate or lorentz iff |aop — Aq is less 
than, equal to or greater than zero. 
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The result follows. □ 

Consider now the case where there are two focal points on each line of S. Thus 
|foP — > 0, and further suppose that these focal points form two continuously 
differentiable surfaces Si and ^2 in M^. Let L be the distance between the focal 
points and if the angle between the normals to Si and ^2 at corresponding points. 

Theorem 5. The distance L and angle ip defined above are given by 

POPO - foCTo I (Top 

Proof. The first of these follows trivially from the fact that the two focal surfaces 
are given by 

9o + (|aop-Ag)^ eo-(|<To|2-A2)i 
y-^ _ u — I ^ — r2 — 



POPO — CQCTO POPO — fOCO 

The line congruence S, by assumption, is not flat, and so we parameterize it by its 
direction ^. To compute the angle if we note that parametric equations for and 

^2 

Z = Zi{^,^) t = ti{^,0 Z = Z2i^,^) t = t2{^,^) 

are obtained by inserting r = ri and r = r2 in equations (1.4). Let vi,i'2 & be 
the directions of the normals to and 52, respectively. Thus, for i = 1,2, 

Vi dzidti — dtidzi 

1 — ViVi dzidzi — dzidzi 

If we introduce, for i = 1, 2, 

a, = dzidti — dtidzi bi = dzidzi — dzidzi 

a straightforward computation shows that 

COS(f = ± j- 

[(6?+4aiai)(6i + 4Q!2a2)]' 

A lengthy computation involving the explicit expressions for aj and bi obtained by 
differentiation of (1.4), yields 

6162 + 2(aia2 + a2ai) = [^dLf - ao{dLf + 2\oidLdL 

+2L{pdL - (3dL) - 4/3^^0 + 4y3Vo - 8Aoi/3y9] 

and 

(6? + Aaiai){bl + ACX2CX2) = -^^—7^ [^o{dLf - cT^{dLf + 2X^idLdL 

+2L{pdL - pBL) - 4/3^^0 + 4^Vo - SXoiPP] ^ 

where we have introduced 

Vko(1 + COV \i^o 
The expression for the angle ip follows. □ 
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6. Reflection off a Cylinder 

Consider a cylinder of radius a, with axis lying along the a;"^— axis in R"^. Then 
we have: 

Proposition 5. The inward pointing normal to such a cylinder is given paramet- 
rically by: 

C = T] = -u e™ 

for {u, v) & R X . The distance of a point p on the surface from the point on the 
normal through p that lies closest to the origin is r = —a. 

Proof. This can be checked by noting that, with the aid of (1.4), the mapping 

(u, v) w), r]{u, v), r{u, v)), with ^, i] and r as stated, yields a parameteriza- 

tion of the cylinder: {u,v) i-^ (— acos w, — asinu , u). Moreover, the oriented normal 
at the point {u, v) on the cylinder is given by the expression in the proposition. □ 

The coffeecup caustic is obtained by finding the focal set of a plane wave reflected 
off the inside of this cylinder. This turns out to be: 

Proposition 6. Consider the reflection off the inside of a cylinder of radius a of 
a line congruence consisting of parallel rays traveling along the x^-axis making an 
angle (3 with the x^-axis. The focal set of the reflected line congruence is a surface 
given parametrically by 



x^ = a cos V I cos^ V — - 



LI \ 



, a ^ 

X = —u cos V cot p, 

2 

forueR and 7r/2 < w < 37r/2. 

Proof. The normal congruence to the plane wave is {^i,r]i), where and jji € C is 
free. Reflecting this off the cylinder given above, wc have by Theorem 3: 

6 = -ae'^" V2 = -^ (ae-- - 2<i - ae^^f ) e^- 

or if we parameterize by the point of reflection /j, = r]o = —ue™ 

We now compute the optical scalars for this line congruence via equations (2.2) and 
flnd that the line congruence is flat (i.e. pp — aa = 0). Thus by Corollary 1 there 
is exactly one focal point on each line, given by, after some computation: 

^ acosf (2cos^ /3 — 1) 

r = — ucosp 

2 sin/? 

Inserting this, along with the expressions for ^2 and r]2 in (1.4) yields the stated 
result. The domain of v must be restricted to half a circle as the incoming rays 
reflect on the inside of only one half of the cylinder. □ 

Wc now consider the focal set formed by reflection of a point source off the inside 
of the cylinder. To this end, the following theorem describes reflection in a cylinder 
as a mapping (^i,r?i)i-^(6,^2): 
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Theorem 6. A ray f^i, rji) intersects a cylinder of radius a lying along the 
iff: 



< a. 



6(1 + 66) 

For such a ray, the reflected ray is 



(6.1) 



6 = -6 



a6(i + 66) 



(«'66(i + 66)' + (6m - 6m)')^ 

J- + 64i 



- 6'7i ± (a'66(i + 66)' + iiivi - ^ivi 



a6(l + 66 



(6.2) 



where exterior and interior reflection are given by the plus and minus signs, respec- 
tively. 

Proof. Consider an incoming ray (^i, r]i). The reflection equations (4.1) and (4.2) 
tell us again that the reflected ray is 



6 = -6e^ 



2iv 



m = -^ (ae-™ - 2<i - ae'''^,f) e^*". 



(6.3) 



(6.4) 



This incoming ray intersects the cylinder iff (c/(4.3)): 

We eliminate u from this equation by combining it with its conjugate and solving 
the resulting equation for v. The solution, which exists iff (6.1) holds, is 

^iv _ 6^1 - 6m ± («'66(i + 66)' + (6^1 - 6m)')^ 

a6(i + 66) 

Substituting this back into the intersection equation we get that 

u = (6m + 6m ± T^(a'66(i + 66)' + (6m - 6m)')^ 

Finally, putting these last two equations into the reflected ray equation (6.3) yields 
the stated result. 

□ 

For multiple reflection we have the following: 
Theorem 7. The If^ reflection of a ray (^i , rji) off the inside of the cylinder is: 



a+i = (-i)'=6 



*ii-(l6l'-*?)^ 



16 



2k 



Vk+1 



i-ll 

6 



where 



6r,i-Mi-|6l')(l6l'-*D^ 
6m - 6m 



161 



2k 



*1 = 



ai(l + |6P) 
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Proof. This follows from iterations of the above theorem once we realize that 
and 

are preserved by reflection in a cylinder. □ 

For a point source at a finite distance, the following theorem describes the focal 
set of the fc*^ reflection: 

Main Theorem 2. 

Consider the fc*'* reflection off the inside of a cylinder lying along 3? -axis with 
radius a of a point source at (—1, 0,0). The focal set of the reflected line congruence 
consists of a surface: 

^ -j^^fc+i^ sinu i + {a^ — Z^sin^ v)i]^''[2klcosvsmv e"*" + (a^ — P sin^ v)i] 

a^'^ ]^kl cos v+ {a^ — P sin^ v) ' ] 

3 _ fc(l - u2)[a2 - /2 _ 2|2sin% + 2klcoav{a^ - P sm^ v)^ 
u\p,kl cos V + (a^ — Psiv? v)^] 

and a curve in the x^x"^ -plane: 

z = {-lf+^ka-^''[lsmv i + ia"- f sm^ v)^f''[l + 2ke'''{a^ - sm^ v)^ 
where z = x^ + ix^, u e K and v is in the domain 

< V < IT for I < a 

and 

— sm^^{a/l) < V < sm~^{a/l) for I > a. 



Proof. Consider a point source lying at {—I, 0, 0) in M^. This line congruence can 
be parameterized by its direction and rji = — Z(l — ^f)/2. The line congruence 
obtained from k reflections of this point source off the inside of a cylinder of radius 
a is given by Theorem 7. 

We then compute the optical scalars of this line congruence parameterized by 
using (2.2). We find, for example, that 

j^^^l k{l + M^) I kl cosv + (a^ — P sin^ v)^] 
'^o 2u(2klcosv + {a^ -Psm^v)^^ 

where = ue'". The similar expression for po shows that the refiected congruence 
is not flat. Thus each line contains exactly two focal points which can be obtained 
by inserting the solutions of the quadratic equation of Theorem 4 into (1.4). The 
results are as stated above. □ 
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7. Discussion 

The focal surface obtained from a plane wave reflected off the inside of a cylinder 
(Proposition 6) is symmetric along the x'^-axis and intersects any plane parallel to 
the x^a;^— plane in a curve. This curve, called a nephroid, is often observed on the 
top of a cup of coffee in the presence of a strong, low and distant light source - 
hence the sobriquet - the coffeecup caustic. Note that this level set is independent 
of the angle /3 of incidence of the incoming light. 

th 

The focal surface of the fc reflection of a point source also has a symmetry: 

Corollary 2. The focal surface generated by the if^ reflection of a point source is 
invariant under translation along the cylindrical axis. 

Proof. This follows from the fact that and x^ of the surface in Main Theorem 2 
are independent of m. □ 



This symmetry is not shared by the reflected wavefront itself - just it's focal sur- 
face. In Figure 7 we illustrate the 1^* reflected wavefront - the lack of translational 
symmetry is obvious. 




FIGURE 7 



In the sequence of pictures in Figure 8 we show the evolution of the level sets of 
the 1^* focal set as the distance of the source decreases. The cylinder is shown by 
the heavy circle. 
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The focal curve lies entirely outside of the cylinder and so is not seen in reality. 
The cross-section of the focal surface inside the cylinder for varying values of I /a 
is illustrated in Figure 9. The parallel wave limit, given by Proposition 6, is also 
indicated with a broken line. 

1 1 1 1 1 




J I I I L 



FIGURE 9 

In Figure 10 we compare the higher reflection caustics for varying values of I /a. 
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Number of reflections 
1 2 3 4 5 




FIGURE 10 

These higher reflection caustics are in fact visible physically. At each reflection 
some of the light intensity is lost and one expects to see a series of overlapping 
caustics of lessening brightness. In fact, the detailed profile of light intensity near a 
caustic varies in ways that geometric optics does not model well. Nonetheless, the 
accompanying plate is a photograph of the caustics formed by a 7cm diameter brass 
cylinder and agrees well with the geometric optics approximation. The photograph, 
which was taken by the first author in collaboration with Grace Weir, shows the 
first and second refiection caustic formed by a light source at l/a — 1 (compare 
with the first two curves on the top row of Figure 10). 
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